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ABSTRACT 

The 4.8 GHz lightcurves of the scintillating intra-day variable quasar J1819+3845 dur- 
ing 2004-5 exhibit sharp structure, down to a time scale of 15 minutes, that was absent 
from lightcurves taken prior to this period and from the 2006 lightcurves. Analysis of 
the lightcurve power spectra show that the variations must be due to the emergence 
of new structure in the source. The power spectra yield a scattering screen distance 
of 3.8 ± 0.3 pc for a best-fit viss = 59 ± 0.5 km s -1 or 2.0 ± 0.3 pc for the scintillation 
velocity reported by Dennett-Thorpe & de Bruyn (2003). The turbulence is required 
to be exceptionally turbulent, with Cff > 0.7 AL^. 1 m -20 / 3 for scattering material of 
thickness AL pc pc along the ray path. The 2004 power spectrum can be explained in 
terms of a double source with a component separation 240 ± 15 fias in 2004. 
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1 INTRODUCTION 

The quasar J1819+3845 exhibits 20-35% rms intensity mod- 
ulations on intra-hour timescales at centimetre wavelengths. 
The variations are due to interstellar scintillation caused by 
plasma thought to be located only z =4-12 pc from Earth 
(Dennett-Thorpe & de Bruyn 2002, 2003, hereafter DB03). 
For z =10 pc the ~30min scintillations observed in the pe- 
riod December to March each year imply an overall source 
size of ~ 60 ^tas, corresponding to 0.39 pc at its redshift of 
0.54 (Ho = 70kms _1 Mpc- 1 , Q. M = 0.27 and fi A = 0.73). 

The source has exhibited variations every time since 
they were first observed in 1999. This is surprising because 
it might be supposed that such a compact source, with 
a brightness temperature exceeding the inverse Compton 
limit, should expand and cease scintillation after several 
months. This problem is common to several persistent ex- 
treme intra-day variable sources including PKS 1519—273 
and PKS 1257-326 (Macquart et al. 2000, Bignall et al. 
2003). At the very least, one would expect such compact 
sources to exhibit some degree of structural variability, made 
evident by changes in the variability characteristics. 

The detection of structural changes in intra-hour vari- 



able sources is complicated, however, by annual cycles in 
their variability time scale (DB03; Bignall et al. 2003). 
Both the amplitude and direction of the scintillation veloc- 
ity change as the Earth's velocity changes relative to the 
scattering medium as it orbits the Sun. This changes the 
scintillation direction, which in turn alters the character 
of the scintillations (i.e. its power spectrum) if the turbu- 
lence in the scattering medium is anisotropic or if the source 
structure is asymmetric. Both asymmetries hinder searches 
for structural variability. These asymmetries are particularly 
important for J1819+3845 because its scintillation pattern 
is anisotropic, with an axial ratio of 15tg* 30 (DB03). 

In this paper we identify the emergence of a new struc- 
ture within J1819+3845 at A6 cm by comparing observations 
taken at identical epochs each year, thus removing ambigu- 
ity introduced by the annual cycle between intrinsic source 
structural evolution and asymmetry in the scintillation pat- 
tern. In §2 we present the lightcurves indicating the emer- 
gence of this new structure. In §3 the lightcurve power spec- 
tra are analysed to deduce its properties. The implications 
of the discovery are summarised in §4. 
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2 OBSERVATIONS 

Figure [1] shows the lightcurves from observations made over 
six epochs in the period 2000-6 at 4.8 GHz with the Wester- 
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and 8. 



Figure 2. Lightcurves of J1819+3845 at 4.8 GHz taken during 
February and March 2004/5. The Feb data have 60s sampling, 
the other two 30 s. 

bork Synthesis Radio Telescope. Details of the data reduc- 
tion procedure are described in DB03. Each observation was 
made at the same time each year so that the scintillation 
velocity vector is identical on all occasions; no lightcurves 
during the period 2001-2 are shown because we do not have 
good quality lightcurves near the pertinent dates. 

The lightcurves exhibit an increase in the degree of 
small timescale structure in 2004-5 compared to those span- 
ning the interval 2000-3. The power spectra of the varia- 
tions develop substantial structure at angular frequencies 
between 0.005 and 0.02 rads , and peak-to-peak varia- 
tions shorter than 15 mins (corresponding to structure on 
a timescale ~ 4 mins in the structure function) are observed 
in the 2004-5 lightcurves. This corresponds to the highly 
significant w ~ 0.005 rads -1 bumps evident in the 2004-5 
power spectra. An additional feature at to ~ 0.008 rads" 1 is 
also evident, albeit at a low power, in the 2005 spectrum. 
Both these bumps are either absent or less evident in the 
2000-3 and 2006 spectra, although several do possess a sig- 
nificant 'knee' near 0.004 rads - . 

Analysis of similar power spectra taken on calibrator 
sources allow us to be confident that the power spectrum 
is not dominated by telescope noise and systematic effects 
down to angular frequencies lo « 0.015 rads -1 (see Mac- 
quart & de Bruyn 2006). Another set of observations taken 
on 26-31 Mar over the period 2001-6 (see Fig.[2]for a com- 
parison between Mar 2004-5 variations with the 21 Feb 2004 
variations) shows a similar recent increase in the amount of 
small scale structure in the lightcurves, but we shall concen- 
trate on the analysis of the February datasets in this paper. 



3 SCINTILLATION MODELLING 

The differences between the character of the lightcurves over 
the interval 2000-6 can be interpreted either in terms of 
(i) changes in the interstellar turbulence responsible for the 
scintillations observed in J1819+3845, or (ii) changes in the 
internal structure of the quasar. This is illustrated by the 
expression for the power spectrum of temporal variations in 
the regime of weak scintillation, applicable to the variations 
of this source at A 6 cm (Codona & Frehlich 1986; Dennett- 



Thorpe & de Bruyn 2000): 




where V(r) is the source visibility measured on a baseline 
r, k — 2n/\ is the wavenumber, viss is the scintillation 
velocity, here oriented along the a;-axis, and 

P pt (q) = 8^r 2 A 2 <J>^ (q) sin 2 , (2) 

is the spatial power spectrum of intensity scintillations that 
would be observed if the source were pointlike. The ampli- 
tude of the line-of-sight-integrated electron density power 
spectrum $jv e , for turbulence located on a thin screen of 
thickness AL at a distance z, is SM = C%AL. The scintil- 
lations are likely to emanate from only a single thin scatter- 
ing screen because only a single screen velocity is required to 
model the annual cycle in the scintillation timescale (DB03) . 

Changes in the turbulence spectrum are disfavoured as 
the cause of changes in the variability characteristics be- 
cause they require a rather contrived modification of the 
power spectrum over only a narrow spatial wavenumber 
range. The appearance of spectral features over the fre- 
quency range u ~ (4 — 10) x 10 -;l rads -1 requires an en- 
hancement of power in the electron density power spectrum 
only over wavelengths 2-kviss /lu ~ (3 — 8) x 10 7 t>5o m (where 
V50 = ?;iss/50kms -1 ). Changes in the variability character- 
istics are more readily interpreted in terms of evolution in 
the micro-arcsecond source structure. Such evolution is ex- 
pected since both the intrinsic flux density and polarization 
are observed to change in this source on > 0.5 yr timescales. 

The ~ 15min variations in the 2004-5 lightcurves re- 
quire the presence of ultra-compact structure within the 
source. However, the detailed shape of their power spec- 
tra depend on both the source structure and the scintilla- 
tion physics. The peaks observed in the 2004-5 power spec- 
tra could be due to either structure in the source visibility 
function or oscillations caused by the sine-squared term in 
cq. ([2]), the 'Fresnel filter'. The latter could become more 
prominent if the source becomes more compact. Oscillations 
due to the Fresnel filter are exacerbated when the turbulence 
is anisotropic, here parameterized by introducing an elonga- 
tion, R, in the electron density spectrum, oriented at an 
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Figure 1. The fias structural evolution of J1819+3845, demonstrated by changes in character of the lightcurves and their associated 
spectra (black points). The red points show the power spectra weighted by a; 8 / 3 (see text for details). The spectra were pre- whitened to 
prevent leakage of low-frequency power associated with the time sampling window function to high frequencies. A noise bias, determined 
from the spectral power near O.lrads -1 , is also subtracted from the spectra (see Macquart & de Bruyn 2006). To reduce the errors, 
at the expense of information at low lj, each spectrum shown is the average obtained by computing separate spectra for each 3-hour 
block of data within each 12-hour lightcurve. The error on each point is half of its mean value and follows a x 2 distribution with 8 
degrees of freedom. It would be appropriate to normalise spectra by the source mean flux density if the entire source were subject to 
scintillation, but as this is not necessarily the case and it is misleading otherwise, the plots here are not so normalised. The 2000-5 
lightcurves were observed with 60 s time resolution while the 2006 lightcurves were observed with 10 s resolution. The source mean flux 
was around 300 mjy over 2003-6 after a fast rise from 125 to 300 mjy between 2000 and 2003. 



angle 9 to viss: 



SM 



(q x cos 9 + q y sin ( 
R 



+ R (q x sin 9 — q y cos 0)" 



Fig. [3] shows the power spectrum for a point source for a 
variety of scattering geometries assuming the index for Kol- 
mogorov turbulence of f) = 11/3. Fresnel oscillations are 
potentially important in the interpretation of J1819+385 
power spectra because its scintillation pattern is anisotropic. 

It is difficult to attribute all the features of the 2004- 
5 power spectra to Fresnel oscillations because the spectral 



peaks change positions between years (see, for example, the 
4mrads _1 feature in the 2004 spectrum). Even within a 
single power spectrum the locations of the peaks are incon- 
sistent with a unique value of tf; numerical integration of 
eq.{T} shows that no combination of rp and 9 reproduces the 
peak locations observed in the 2004-5 spectra. Nonetheless, 
detailed fitting described below suggests that the first peak, 
at 1.7mrads , exhibited by many of the spectra, can in- 
deed be attributed to the first peak of the Fresnel filter, and 
implies a screen distance ~ 5 v$o pc. 
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The inability of the Fresnel filter to account for the 
peaks in the 2004-5 spectra leads us to consider the addi- 
tional influence of source structure. We first discuss sim- 
ple arguments to deduce general source properties from the 
power spectra before reproducing them quantitatively using 
a specific source model. 

We examine the size of the structure implied by the 
fast variations evident in the 2004-5 lightcurves. The scintil- 
lation power expected for a point source declines as lu~ 8 ^ 3 
for uj > 2.5 x 10" 3 rads _1 and (3 = 11/3 (viz. eq. Q and 
Fig-® - F° r a point-like source the visibility, V(r), and hence 
w 8 / 3 P(uj), would be flat. Consider the u 8,/3 -weighted power 
spectra shown in red in Fig.[T] Fluctuations at uj depend on 
\V(r Pa ijJTp/viss)\ so the frequency at which the red points 
depart from a flat line marks the point at which the source 
structure dominates the shape of the power spectrum. This 
yields a rough estimate of the source size. A sharp decline in 
u; 8//3 P(oj) at a scale too indicates the presence of structure 
down to an approximate scale 

q . g 3 f uo \~ 1 ( VISS \ ( z \ 1 as (3) 
VO.Olrads^ 1 / V 50km s" 1 / \ 10pc/ MaS ' 

For instance, the 0.01 rad s _1 break in the 26 Feb 2006 power 
indicates structure down to ~ 7t>so z^ 1 fias if the index of 
the turbulence spectrum is j3 — 11/3. If the turbulence spec- 
trum were shallower there is less source structure at high 
angular resolution, and if steeper there is more. 

Since the bumps in the power spectra cannot be repro- 
duced in terms of Fresnel oscillations, we consider the sim- 
plest source model capable of introducing additional spectral 
peaks. A double source, with component flux densities 7i 
and I2 and sizes 81 and 82, with brightness power spectrum 

nrr \|2 t 2 -k 2 r 2 B 2 , T 2 -k 2 r 2 Bi 

\V(r)\ = I x e 1 + I 2 e 2 

+2I 1 l2e- k2r2(e i +02)/2 cos(£:A6> ■ r), (4) 

induces oscillations in the scintillation power spectrum 
whose spacing depends on the component separation, Ad = 
A8 (cos q, sin a), where a is the angle the vector makes 
with viss- The oscillations are partially damped if the mean 
square of the component angular sizes exceeds ~ A8. 

We performed fits to two power spectra up to uj = 
0.015 rad s _1 , the point beyond which the power dips be- 
low the level at which we are able to confidently measure 
flux density variations. The 26 Feb 2006 spectrum, closely 
resembling the 2000-3 spectra, was used as representative of 
spectra with the simpler structure, while the 21 Feb 2004 
spectrum is representative of data in which additional fast, 
small amplitude variations are prominent. The fit was per- 
formed to both spectra simultaneously to ensure the scatter- 
ing parameters, SM , z, R, Diss and 8, were common to both. 

The best-fitting parameters are listed in TableQ] Dou- 
ble component models were fit to both spectra. Fig. [4] shows 
that the models reproduce the essential characteristics of 
the power spectra well. Two models, Ml and M2, were fit- 
ted for. They differ in that Diss and a were fixed in M2 to 
values suggested by the annual cycle in the source's variabil- 
ity timescale (DB03). In fixing a instead of 8 we assume that 
the source structure, rather than anisotropy in the turbulent 
scattering plasma, dominates the anisotropy of the source's 
scintillation pattern. This is justified by the low value of 
R derived in all fits (see below). However, we caution that 



the large anisotropy of the scintillation pattern renders the 
solutions for diss and a degenerate in the annual cycle fit 
described in DB03. Thus it is not clear that the "best-fit" 
parameters suggested in DB03 do indeed represent the cor- 
rect solution. 

The quality of the fit was degraded by fitting the spec- 
tra with a single source component. However, we note that 
the best-fit flux density of the second component is only 
2-3<T from zero. A fit without the second component in the 
2006 data yields a reduced \ 2 of 1.2. The double-component 
model considered here is nonetheless strongly justified by 
polarization observations: time delays are observed between 
the polarized and unpolarized variations of this source over 
the entire period 2000-6, demonstrating that the source is 
always comprised of at least two components, one of which 
is likely highly polarized (Macquart, de Bruyn & Dennett- 
Thorpe 2003). There was no improvement in the fit qual- 
ity by allowing (5 to vary from its assumed value of 11/3; 
the reduced \ 2 f° r the fit was 1.02 and yielded /3 = 3.87. 
However, other possible islands of low x 2 conceivably exist 
besides those detected by our fitting program and may yield 
acceptable fits with values of f3 ~ 5. 

The R — 2.9 ± 1.0 turbulence anisotropy derived from 
the fit is smaller than the 15tg* 30 anisotropy of the scin- 
tillation pattern itself (DB03), but is not inconsistent be- 
cause the latter is also influenced by asymmetry in the 
source structure. The value of R is particularly sensitive to 
the low-frequency portion of the power spectrum. Two low- 
frequency points, obtained by recomputing the power spec- 
trum from the lightcurves subdivided into 6 hr blocks, were 
therefore added to the spectra employed in the fitting (see 
FigHJ). This additional information is obtained at the ex- 
pense of poorer error estimates associated with these points. 
The derived anisotropy ratio should therefore be regarded 
with some caution; we expect to derive a better estimate in 
future from longer duration observations. 

There is a degeneracy between SM and the compo- 
nent flux densities; only the product SM I 2 affects the spec- 
tral amplitude. Since the sum of the component flux den- 
sities cannot exceed the 281 and 295 mjy mean flux densi- 
ties observed in the 2004 and 2006 lightcurves respectively, 
the highest possible values of 7i and I2 in the 2006 spec- 
trum in Ml are 236 mjy and 60mJy respectively, requiring 
SM > 2.1 x 10 16 m" 17/3 . This places a lower bound on the 
turbulent amplitude of C% > 0.68 (AL/1 pc)" 1 m" 20/3 . For 
comparison, Macquart & de Bruyn (2006) deduced SM = 
(1 - 1.4) x 10 17 m~ 17/3 for z = 4pc in a fit to the refractive 
scintillation at 1.4 GHz. A Virginia Tech Spectral-Line sur- 
vey measurement (see |http://www.phys.vt.edu/~ha lpha/ ) 
of the Ha intensity at the position of J1819+3845 implies 
SM = 1.1 x 10 17 T 4 a9 /~ 2/3 e 2 /(l + e 2 ) m" 17 / 3 , where the gas 
temperature is T = 10 4 T4 K, lo is the turbulence outer scale 
in AU and e 2 = {(5n e ) 2 ) jrl 2 is the normalised electron den- 
sity variance. This suggests lo < 9 AU if e ~ 1 and T4 < 10. 

A weaker degeneracy between the screen distance and 
component angular sizes allows decreases in z accompanied 
by corresponding increases in all measured angular sizes to 
yield fits of similar quality. However, changes in z do ad- 
versely affect the fit because they also alter the positions of 
the Fresnel oscillations. 
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Figure 4. A fit to the 21 Feb 2004 and 26 Feb 2006 power spectra 
with the parameters of Table 1 (Ml blue line, M2 green line). The 
red line represents Ml but with R = 3.5. The higher anisotropy 
fit provides a better fit to the low-frequency part of the 2004 
spectrum, but a worse fit overall with a reduced \ 2 of 1.3 vs. 1.0. 



Ml 



reduced x 2 
Fit Parameter 



39.1/39d.o.f. 



1.0 



M2 
45.8/42d.o.f. 



1.1 



SM(10 17 m 


-5.67) 


2.6 ±0.2 


2.4 ±0.2 




z(pc) 




3.8 ± 0.3 


2.0 ±0.3 




viss (km s~ 




59.0 ±0.5 


42.5* 




R 




2.9 ± 1.0 


1.9 ±0.6 




9 (rad) 




2.26 ±0.02 


0.9 ±0.2 






2004 


2006 




Ml 


M2 


Ml 


M2 


h (mjy) 


42 ± 7 


69 ±30 


67 ± 11 


100 ± 45 


h (mjy) 


11 ±4 


9±4 


17 ± 12 


25 ± 13 


A6 (/ias) 


240 ± 15 


220 ± 20 


565 ± 15 


556 ± 16 


a (rad) 


0.00 ±0.03 


1.07* 


2.4 ±0.3 


1.07* 


8i (fi&,s) 


16 ±1 


18 ±2 


20 ±3 


26 ±4 


2 (pas) 


16 ±2 


0±500 


7±4 


9 ± 3 



Table 1. Fit parameters with 1-tr errors for the 2004 and 2006 
power spectra derived using the CERN MINUIT minimization 
package, which estimates errors based on the second derivatives 
of x' 2 ■ Fit parameters are defined in the text. Fitting models Ml 
and M2 differ only in that parameters marked by asterisks are 
held fixed in M2. The redshift-corrected brightness temperature 
of component 1 in Ml is 1 X 10 13 K. 



4 CONCLUSION 

During 2004-5 fast, ~ 15 min, variations emerged over the 
~ 30 min variations normally observed in the lightcurve of 
the scintillating quasar J1819±3845. In 2006 the source re- 
verted to variations similar to those in 2000-3. The changes 
are best explained in terms of evolution in the source struc- 
ture. A double-component source is capable of explaining 
the power spectra during 2004-5. The spectral peak at 
w~2x 10~ 3 rads -1 corresponding to the large amplitude, 



slow variations and visible in all spectra, is explained by 
oscillations of the 'Fresnel filter', from which a 5vsopc dis- 
tance to the scattering medium is deduced. Holding the scin- 
tillation velocity fixed at the value derived by DB03 yields 
z = 2.0T0.3, while fitting also for wiss yields 3.8T0.3 pc. The 
fast variations in 2004-5 are best attributed to oscillations 
in the source visibility function caused by a double-source 
structure with a component separation of 240 ± 15/xas in 
2004. From a fit to the 2006 spectrum one infers a sepa- 
ration of 565 ± 15 /xas. This implies an apparent expansion 
speed of 3.4 ± 0.3 c over the two years separating the two 
observations. However, we caution that it is also possible to 
fit the 2006 spectrum with in a single component, in which 
case no expansion speed can be derived. 

The emergence of these fast variations preceded other 
changes in the source, namely an increase in the modulation 
index and then the intrinsic flux density. The appearance of 
new structure may be connected with the diffractive scin- 
tillation reported at 21cm (Macquart & de Bruyn 2006), 
but as absolute astrometry is not possible using scintillation 
techniques one cannot determine the angular offset between 
the components responsible for 6 and 21 cm variations. 

The proximity of the scattering screen and the large 
amplitude of the intensity variations requires the scattering 
turbulence to possess an amplitude C% > 1.7m~ 20 ' 3 if lo- 
calized to a region of thickness no more than 0.4 pc, 10% 
of its distance. This exceeds C% values deduced from most 
scattered pulsars (Cordes et al. 1988) by over two orders of 
magnitude. 
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